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Dynamical properties of an impurity spin coupled symmetrically to sublattices of ordered 2D 
Heisenberg quantum antiferromagnet (i.e., frustrated impurity spin) are discussed at T > (exis- 
tence of a small interaction stabilizing the long range order at T 7^ is implied). We continue our 
study on this subject started in Phys. Rev. B 72, 174419 (2005), where spin- 7^ defect is discussed 
and the host spins fluctuations are considered within the spin- wave approximation (SWA). In the 
present paper we i) go beyond SWA and ii) study impurities with spins S > 1/2. It is demonstrated 
that in contrast to defects coupled to sublattices asymmetrically longitudinal host spins fluctuations 
play important role in the frustrated impurity dynamics. We show that the effect of the host system 
on the defect is completely described by the spectral function as it was within SWA. The spectral 
function, that is proportional to to 2 within SWA, acquires new terms proportional to ui 2 and ujT 2 
originating from longitudinal host spins fluctuations. It is observed that the spin- i impurity suscep- 
tibility has the same structure as that obtained within SWA: the Lorenz peak and the non-resonant 
term. The difference is that the width of the peak becomes larger being proportional to / 2 (T/J) 3 
rather than / 4 (T/J) 3 , where / is the dimensionless coupling parameter. We show that transverse 
static susceptibility acquires a new negative logarithmic contribution. In accordance with previous 
works we find that host spins fluctuations lead to an effective one-ion anisotropy on the impurity 
site. Then defects with S > 1/2 appears to be split. We observe strong reduction of the value of the 
splitting due to longitudinal host spins fluctuations. We demonstrate that the dynamical impurity 
susceptibility contains 2S Lorenz peaks corresponding to transitions between the levels, and the 
non-resonant term. The influence of finite concentration of the defects n on the low-temperature 
properties of antiferromagnet is also investigated. Strong spin-wave damping proportional to n/ 4 |w| 
' O ■ is obtained originating from the non-resonant terms of the susceptibilities. 
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I. INTRODUCTION 



The problem of impurity spins in 2D Heisenberg quantum antiferromagnet (AF) has been attracted much attention 
in last two decades ,L2iiMiiii£iL2i2ii2iiiiI2iiiiI4 The majority of works are devoted to defects coupled asymmetrically to 
sublattices of AF. Among such impurities are added spin coupled to one host spin, substitutional spin and vacancy 
t-h ' that is the particular case of the added spin with the coupling strength g — > 00. Such works are stimulated by a 
\& I variety of experiments studying cuprates with magnetic and nonmagnetic impurities. Defects coupled symmetrically 
to sublattices of AF (see Fig.^l are much less studied because such objects are much rare in occurrence. At the same 
time they have different propertiesii*^ and are also of interest. For instance, attempts were made to model holes in 
CuO-planes of some high-T c compounds being in antiferromagnetic phase by spin-i impurity coupled symmetrically 
to two neighboring host spins 

In our recent paper— (hereafter referred to as I) we studied dynamical properties of impurity coupled symmetrically 
to two neighboring spins in 2D AFs at T > 0. A technique was proposed based on Abrikosov's pseudofermion 
technique^ that allows to discuss dynamics of an impurity with arbitrary spin value S. Meanwhile we focus in I on 
O ■ the particular case of spin-i defect. Our aim was to find the impurity dynamical susceptibility xi 1 - )- We assumed 
that the impurity dynamics is governed by interaction with spin waves. Then, two kinds of the host systems are 
considered in I: i) ordered 2D AF in which the long range order at T ^ is stabilized by a small interaction (for 
definiteness interplane interaction) 77 <C J, where J is the coupling constant between the host spins; ii) isotropic 2D 
AF. We bared in mind that in isotropic 2D AF only spin waves are well defined with energies much larger than Jsa/£, 
where a is the lattice constant, s is the value of host spin, and £ oc exp(const/T) is the correlation length. It was 
obtained within the spin-wave approximation (SWA) that, similar to the spin-boson model, the effect of the host 
system on the defect is completely described by the spectral function which was assumed to be proportional to ui 2 in 
calculations of x(w). We demonstrated that the spectral function is proportional to u 2 in 2D AF and in the ordered 
quasi-2D AF at u> ^> Jsa/£ and uj ^S> 77, respectively. Notice that only transverse host spins fluctuations contribute 
to the spectral function within SWA. Then, only processes of absorption/emission of one spin wave by the impurity 
is taken into account within SWA. 

The calculations of xi^) were performed in I within the fourth order of the dimensionless coupling parameter / oc 
g/J. We shown that the transverse impurity susceptibility x_i_(w) has a Lorenz peak with the width r oc f A J{T/jf 
that disappears at T = 0, and a non-resonant term. The imaginary part of the non-resonant term is a constant 
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independent of T at \u>\ 3> r and the real part has a logarithmic divergence of the form / 2 ln([w 2 + r 2 ]/J 2 ). The 
longitudinal susceptibility has only the non-resonant term which differs from that of X±( UJ ) by a constant. 

The fact that the spectral function in 2D AF is proportional to oj 2 only at oj ^> {77 or Jsa/£} leads to the following 
restriction on the range of validity of the results obtained: max{_T, \oj\} ^> {77 or Jsa/£}. Notice that in the case of 
isotropic 2D AF at T ^ this relation determined the range of validity of our assumption that the impurity dynamics 
is governed by interaction with spin waves. 

It was observed that the static susceptibility x(0) has the free-spin- like term 1/(4T) and a correction proportional 
to f 2 ln( J/T). Quite a different behavior of x(0) was obtained in the regime of T <C \g\ for asymmetrically cou- 
pled defects^* 6 - the classical-like term S 2 /(3T) and the logarithmic correction proportional to ln(J/T) rather than 
PHJ/T). 

The influence of the finite concentration of the defects n on the low-temperature properties of 2D AF was also 
studied in I. It was shown that correction to the square of the spin-wave spectrum is proportional to nf 2 x±( u >)- Then 
two regimes were considered, \u>\ > wo « f 2 T 2 /J and \u\ <C ll>o, at which, respectively, non-resonant and resonant 
parts of Xi-( w ) prevail. In the most interesting case of \lo\ 3> wo we found the logarithmic correction to the spin-wave 
velocity and an anomalous damping of the spin- waves proportional to n/ 4 |w|. Such damping was obtained also in 
Refsii& where asymmetrically coupled impurities were studied in 2D AF. 

In the present paper we continue our discussion of symmetrically coupled defects in 2D AFs and i) go beyond SWA 
and ii) study defects with 5* > 1/2. In particular, longitudinal host spins fluctuations come into play outside SWA. 
They lead to two kinds of processes: absorption/emission of two magnons by the impurity and scattering of one 
spin wave on the impurity. We find that in the case of isotropic 2D AF at T ^ our assumption that the impurity 
dynamics is governed by interaction with spin waves is wrong at any T and to. This circumstance manifests itself in 
the fact that spin waves with wavelength of the order of £ are important in the processes of absorption/emission of 
two magnons and scattering of one magnon. Thus, our approach can be applied to ordered 2D AFs only. 

It is demonstrated that longitudinal host spins fluctuations lead to two important contributions to the spectral 
function which should be taken into account. One of them is proportional to T 2 to/s and corresponds to scattering of 
one spin wave on the impurity. Another contribution has the form uj 2 v(T)/ s, where v(T) is a series in (T/s) ln(T '/ [srj\) , 
and originates from two processes: emission or absorption of two spin waves by the impurity and the scattering of 
one magnon on the impurity. 

It is demonstrated that spin-i impurity susceptibility Xx(w) has the same structure as that obtained within SWA: 
the Lorenz peak and the non-resonant term. The difference is that the width of the peak r acquires large correction 
of the first order of f 2 : r oc / 2 J(T/J) 3 . This term dominates if the following condition is fulfilled that does not 
depend on s: \g\/ J <C 5. Besides, some constants in the expression for Xx(^) acquire thermal corrections. The width 
of Lorenz peak in XH^) is zero within our precision. We observe that longitudinal host spins fluctuations lead to a 
new logarithmic correction to XJ-(O) proportional to — (f 2 /[Js]) \n(T/[srj\). 

Corrections to the spin-wave spectrum have the same form as those obtained within SWA. The difference is that 
some constants acquire thermal corrections, r has another form and the quantity ojq separating two regimes becomes 
larger being proportional to T 2 / J rather than f 2 T 2 /J. 

It should be pointed out that the obtained remarkable influence of longitudinal host spins fluctuation on frustrated 
impurity spin dynamics is quite unusual. There is no such influence in the case of asymmetrically coupled defects. 
We show below that longitudinal fluctuations give a negligibly small correction to the spectral function and there 
is no reason to expect any perceptible influence from them. Thus, we confirm that the frequently used T-matrix 
approach 1 - 16,17 basing on Dyson-Maleev transformation and dealing with only bilinear part of the Hamiltonian is 
appropriate for discussing the asymmetrically coupled defects. 

We find also x( w ) f° r symmetrically coupled impurities with S > 1/2. It is well-known that host spins fluctuations 
lead to an effective one-ion anisotropy on the impurity site in 3D AF: —C(T)S 2 g 2 / J, where C(T) > Oi£il& Then 
defects with S > 1/2 appear to be split and magnetization of sublattices is a hard axis for them. We obtain the same 
one-ion anisotropy for defects in 2D AF. In particular we observe large reduction of the value of C(T) obtained within 
SWA by 1/s-corrections stemming from longitudinal host spins fluctuations. We show that thermal corrections reduce 
C{T) further and can change the sign of C(T) at T < T N . 

Then we find that if T is greater than the value of the defect splitting, transverse dynamical susceptibility contains 
2S Lorenz peaks corresponding to transitions between impurity levels and a nonresonant contribution. Widths of the 
peaks are proportional to f 2 J(T/J) 3 as in the case of S = 1/2. When T is lower than the distance between nearest 
impurity levels (~ C(T)g 2 / J) there are only peaks corresponding to transitions between low-lying levels. Longitudinal 
susceptibility has only the nonresonant term. At T ^ C(T)g 2 S 2 / J static susceptibility has the same structure as 
that of S = 1/2: term S(S+ 1)/(3T) and a logarithmic correction. If T < C(T)g 2 /J, x(0) of integer and half-integer 
spins behave differently. The imaginary part of the nonresonant term is a constant that leads to abnormal spin-wave 
damping proportional to n/ 4 |w|, as in the case of spin-i impurity. 
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It is explained in I that in the case of finite concentration of defects the results obtained within our method have 
a limited range of validity. This is because the interaction with defects modifies host spins Green's functions and the 
bare spectral function acquires large corrections at energies around those of the Lorenz peaks. The problem should 
be solved self-consistently to obtain correct expressions for the impurity susceptibility and the spin-wave spectrum 
around the resonance energies. Corresponding consideration is out of the scope of the present paper. This large 
renormalization of the spectral function can be interpreted as a resonance scattering of spin waves on impurities. 
Similar situation exists in the case of asymmetrically coupled defects (see discussion in Sec. 0). As a result in the 
vicinity of the resonance energies one might expect strong deviation from linearity of the real part of the spin-wave 
spectrum and an increase of the spin- wave damping. 

The rest of this paper is organized as follows. The model and diagrammatic technique are discussed in Sec. [n] 
Pseudofermion Green's functions and pseudofermion vertex are calculated in Sec. IIIII The impurity dynamical sus- 
ceptibility is derived in Sec. lIVI Influence of the defects on the spin- wave spectrum is considered in Sec.[V] Section IVll 
contains our conclusions. There are two appendixes with details of calculations. 



II. TECHNIQUE AND GENERAL EXPRESSIONS 



A. Technique 

We discuss Heisenberg AF with impurity spin S coupled to two neighboring host spins si and S2 describing by the 
Hamiltonian 

n = j^s iSi + 0S(si + S2 ). (i) 

It should be stressed that one can consider other symmetrically coupled defects (e.g., those coupled to four host 
spins) on the equal footing. The results would differ by some constants only. It is convenient to use Abrikosov's 
pseudofermion technique^ and to represent the impurity spin as S = ^2 nm &„S nm 6 m , where n and m are the spin 
projections, 6^ and b n are operators of creation and annihilation of pseudofermions. We calculate below impurity 
susceptibility x(o>) using diagrammatic technique. First diagrams for x( w ) an d a graphical representation of the result 
of all diagrams summation are shown in Fig. [21 Thin lines with arrows in the picture represent the bare pseudofermion 
Green's functions: G^ m ,(iu> n ) = 5 mm >(iu) n — A) -1 , where A is the chemical potential of pseudofermions that should 
be tended to infinity in the resultant expressions. Wavy lines in the picture denote magnons Green's functions. As 
usual, diagrams with only one pseudofermion loop should be taken into account because each loop is proportional to 
the small factor of e _A / T . 

Within SWA discussed in I only one virtual magnon can be emitted or absorbed in each vertex. Beyond SWA 
one has to consider diagrams in which some vortexes contain two and three magnons lines. Then, within SWA 
products containing many operators of host spins reduce to products of two-spins Green's functions while there is no 
such simplification beyond SWA. Meanwhile we find below that within the first order of g 2 , the order in which all 
calculations of the present paper are done, diagrams with single and double wavy lines that do not contain vertexes 
with three wavy lines give the largest contributions. These are diagrams shown in Figs. |21 and [3J Fig. Ogives examples 
of diagrams that should be discarded. A double wavy line corresponds to longitudinal susceptibility of the host spins. 
Thus, the diagrammatic technique used in I can be applied for the present task with minor modification: single wavy 
line and double wavy line correspond to 7r _1 A r (a;)ImA /J1/ (ix') with //, v — x,y and it N(u>)JxaA. zz (u>), respectively, 
where N(u>) — (e"/ T — I)" 1 is the Planck's function, 

/>oo 

A^(u,) = -i dte*-*([<(t) + ^(t),^(0)+ fl S(0)]) (2) 



and (...) denotes the thermal average; frequencies of functions ImA /J!y (o;) should be taken so as they are contained 
in arguments of G^-functions with positive sign; integration over all frequencies of functions ImA M „(w) is taken in 
the interval (—00,00). ImA M „(w) is referred to throughout this paper as spectral function. 

It is shown in Appendix 1X1 that the imaginary part of A m1/ (cj) given by Eq. has the following form in the ordered 
2D AF at \u\ > rj: 

ImA^H = -^(^) 2 Bgn(w)A(w)d^ - (§) A(w)^^ Z) (3) 

= (5^(1 - 5 ^ H —SpgSvg, (4) 
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where A and B are positive constants given by Eqs. i|A4(l and i|A8|) . respectively, that are independent of s and which 
dimensionality is inverse energy, is a characteristic energy for which we have Eq. I|A3[) and A(w) is a cut-off function 
that is equal to unity at \w\ < 9 and it drops rapidly to zero outside this interval. Within the first order of 1/s we 
have 

r(T) = l, <T) = ^£^), (5) 

where A*" is the number of spins in the lattice and ek is the spin-wave energy which is equal to \/8sJk at small k. The 
logarithmic infra-red singularity in the expression for v(T) is screened by the interaction stabilizing the long range 
order at T ^ 0: v(T) = T/ (ttsJ) ln(T/ (si])). Notice that only the term proportional to Aw 2 d M1/ (l — 5^ z ) was obtained 
in I within SWA. It originates from transverse host spins fluctuations and corresponds to emission or absorption of one 
spin wave by the impurity. The remaining two contributions to ImA Al ^(ti>) in Eq. Q are of the next order of 1/s. They 
stem from longitudinal host spins fluctuations. The last term in Eq. is proportional to T 2 lu and corresponds to 
scattering of one spin wave on the impurity. This term can be larger than u> 2 . The correction proportional to v(T)uj 2 
originates from two kinds of processes: emission or absorption of two spin waves by the impurity and scattering of 
one magnon on the impurity. To determine range of values of v(T) we remind the well-known formula for the average 
z component of the host spin 

V ' N ^ 2e k N ^ e k y > 

The second and the third terms in Eq. © are equal approximately to 0.2 and T/(2irsJ) \ii(T/[st]]), respectively. The 
third term must be much smaller than s within spin-wave approach. Comparing it with v(T) given by Eq. (0 one 
concludes that v(T)/(2s) is much smaller than unity. 

As it is explained in Appendix lAl higher order 1/s-corrections give contributions proportional to products of w 2 
(ivT 2 ) and powers of [T/(sQ)] ln(T/[srj\). Hence, r(T) and v{T) appear to be series in powers of [T/(sQ)] ln(T/[srj\). 
We restrict ourself in this paper by the first 1/s-correction to r(T) and v(T). 

Term in Eq. J2J proportional to v(T)uj 2 will result in the temperature corrections to some constants in x( w ) obtained 
within SWA. It will give also new logarithmic contribution to x±(0). The last term in Eq. © will lead to the large 
renormalization of the Lorenz peaks widths. 

We show in Appendix 1X1 that one can lead to expressions (Oh® considering isotropic 2D AF at T ^ 0. However 
there is no parameter in this case to screen infra-red singularity of the function v(T). This singularity signifies 
that host excitations with wavelengths greater than the correlation length £ play an important part in the impurity 
dynamics. Then spin- wave formalism is inadequate in this case and our results are applicable for ordered 2D AFs 
only. 

B. Dynamical susceptibility of the impurity 

We have for the dynamical susceptibility after the analytical continuation from the discrete frequencies to the real 

14.20.21.22 



— A/T 

XP (u) = — — - I dxe- x/T Ti{P[G(x + uj)T+ + (x + LJ,x)G(x)-G*(x)r p -(x,x--uj)G*{x-uj) 



axis 

- G(x + uj)r+-(x + uj,x)G*(x) +G(x)T+-(x,x- w)G*(x-w)]}, (7) 

where A/" is the number of pseudofermions that is proportional to e _A / T , P is a projection of impurity spin, G(w) is 
the retarded Green's function, the trace is taken over projections of the impurity spin and signs at superscript of Tp 
denote those of imaginary parts of the corresponding arguments [e.g., Tp~(x,y) = Tp(x + iS, y — ifi)]. An energy shift 
by A has been performed during the derivation of Eq. (7J. As a result the Fermi function (e^ x+x ^/ T + has been 
replaced by e~( x+x ^ T and the functions G and Tp no longer depend on A. These are those functions we calculate 
in the next section by the diagrammatic technique. It is clear that the bare pseudofermion Green's functions in this 
case are G^ m ,(u) = 5 mm >/bj. 

We derive below analytical expressions for the dynamical susceptibility of the impurity. Perturbation theory is used 
for this purpose according to the interaction. It can be done if the dimensionless constants 

2 g 2 A 2 g 2 Br(T) 

f = _ and h = (8) 
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are small. Two constants / and h both proportional to \g\ are introduced to distinguish contributions from the terms 
in Eq. J3J proportional to w 2 and T 2 uj. Notice that h = within SWA. Using Eqs. (|A3() . (|A4I) and l|A8|l we have 
f 2 = (g/ J) 2 y^/ (4s) and /i 2 = (g/ J) 2 it 2 /s 2 . Then g can be even greater than J at large enough s. 



III. PSEUDOFERMION GREEN'S FUNCTION AND THE VERTEX 



A. Green's function 



We turn to the calculation of the pseudofermion Green's function G mn {u>). The Dyson equation for it has the 
following form: LoG mn (uj) — S mn + ^mq{w)Gq n {uS), where E mg (w) are matrix elements of the self-energy. It is easy 
to show that matrices S(w) and G(w) are diagonal: T, mn (uj) = 5 mn Yi n (uj) , G mn (uj) = 8 mn G n (uj) = S rnn /[uj — S n (w)]. 
As we demonstrate in I, there is further simplification in the case of two-level impurity: £(w) and G(w) are proportional 
to the unitary matrix. 

Diagrams of the order of g 2 for £ n (w) are shown in Fig. [3] Let us represent the Green's function in the form 

G„M = - l - Z ^\ (9) 

U) + C n + l"f n (L) + C n ) 

where Z n (uj) and 7n(w) are some functions, 7n(w) are real ones, and c n are constants. Using Eq. we have in the 
first order of g 2 

z m M _ [°° 4, US> (10) 

7rB J_ 00 X + U) + C n + l-f n {X + OJ + C n ) 

7 «H = ^ 2 4 2 ^|^[i + iv(c)]AH, (ii) 

c « = g 2 ^ 5+i )- 2 r te&^+t* r dx m^M (12) 

71" '-oo ^ ^ '-oo X 

where A±(lo) — A xx (lu) = Ay y ((v), the principal value of the integrals is implied in Eq. 112|) and 

4 1 ] = {S»S") nn d nv = S(S + 1) - (l - ^pj n 2 , (13) 

4 2 i = (S z S z ) nn = n 2 . (14) 

Summation over repeated Greek indices is understood in Eq. I|13|) and below. The logarithmic divergence in expression 
(II Ul) at real to is screened by the term i"/ n (x + lu + c„ ) in the denominator. We find that higher order diagrams both in 
g and 1/s (e.g., that shown in Fig.0](a)) give negligibly small contributions to the self-energy part. It was obtained 
in I that within SWA the first nonzero contribution to 7 n (^) is of the order of f 4 cx g A . As is seen from Eq. Ijllfl . 
longitudinal host spins fluctuations lead to the correction of the order of h 2 oc g 2 . It is important to note that 7 n (^) 
is the constant at M <C T: 



7n ( W - 0) « r Qn = fc 2 Jjg>0 . (15) 

Notice that within SWA this constant shows the same T-dependence^ r^ wa) oc f 4 T 3 /0 2 . 

Using Eqs. (|§|l- H12|) the number of pseudofermions M can be calculated up to the first order of g 2 with the result 

1 t°° 

N = -J2~ dxN(x)lmG n (u) « e~ A / T ^ e c " /T . (16) 

Notice that terms proportional to g 2 cancel each other in the right part of Eq. (|16fl . 

It should be stressed that one can not use in Eq. (fT~2f) expression © for imaginary part of A uv (x) because large x are 
essential in the integrals. We have to use exact expressions for ImA Mt/ (x) to find constants c: Eq. I|A2(1 for ImAi(j:), 
Eq. (|A5|) for 1/s-correction to ImAi(i), and Eq. (|A7|) for ImA Z2 (x). As a result one obtains up to unimportant 
constant independent of n 

cl^-(c-^^)4 = -C(T), 2 ^, (17) 



T 
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where 



J r 1 +cos(kR 1 2 ) 
(2tt) 2 J Jo + J k 



C = T^rJ r , ' « 0.25, (18) 



(2tt) 2 J ~"sJ 



C± = C ( 1 - y— ^ / dk-f- ) « 0.04, (19) 



, m sJ f ji. l-cos([k 1 +k 2 ]R 12 )( S Jo) 2 + S 2 J kl J k2 -e kl 6 k2 , r /4sJ f „N(e u ) 



(27r) 4 7 " Ckifka Cki+Cks V( 27r ) 2 J £ k 

(2tt) 2 J e k 



0.lB + £i[*^, (20) 



where R12 is a vector connected host spins coupled to the impurity, J k = 2J(cosk x + cosk y ), e k = s^f Jq — J k is 
the spin-wave energy and integrals are over the chemical Brillouin zone. The constant C is of the zeroth order of 
1/s. It originates from transverse host spins fluctuations. Constants C± and Cy(T) are first 1/s-corrections stemming 
from transverse and longitudinal fluctuations, respectively. It is seen that their sum is of the order of C at small s 
and it should be taken into account. Temperature corrections are considered in C\\(T) only because those to C± are 
much smaller. We note that the thermal correction to Cy(T) coincide with that to (s z ) given by Eq. ©. It should be 
pointed out that this correction reduces the value of C(T) and can change its sign at T < T/y. 36 

Let us discuss the physical meaning of constants c given by Eq. (fTTj) . They describe splitting of the impurity caused 
by the host spins fluctuations. The distances between the levels are given by differences between corresponding 
constants c. One infers from Eq. I|16|) that at T much smaller than the distance between the nearest levels only 
those with maximum c are populated. As is clear from Eq. I|17|l . these are levels with the smallest in absolute 
value projections on the quantized axis. Then, the impurity splitting can be described by the effective Hamiltonian 
—S 2 C(T)g 2 1 J that is the effective one-ion anisotropy. Such anisotropy was obtained within SWA for frustrated defects 
in 3D AFiSiia 

It is also seen from Eq. i|17fl that spin-i defect remains degenerate because ci = c_i. We have taken advantage of 
this circumstance in I and discarded constants c attributing them to renormalization of the chemical potential A. 



B. Pseudofermion vertex 



Let us turn to the consideration of the pseudofermion vertex Tp(x + lu,x). First diagrams for this quantity are 
presented in Fig. It is shown in I that Tp mm i (x + ui, x) is proportional to P mm > for 5 = 1/2. However in general it 
has a more complicated matrix structure. 

As is seen from Eq. (JTJl, we need four different branches of Tp(x + cu, x). It is clear that T ++ = (r - ~)* and within 
the first order of g 2 one has 



j-2 roo 

r £L(* +«,!) = Pmn + ^ E SU miPi q S» n d^ / dyy\y\N(y)A(y)G l (x + y + cj)G q (x + y). 

l,q J -°° 



(21) 



It is seen from Eq. I|21ll that the poles of G-functions are on the one hand from the real axis. Hence, the second term 
in Eq. 1)21(1 is small compared to the first one and we can restrict ourselves by this precision. The term proportional 
to h 2 is discarded in Eq. I|21|l being negligibly small due to the factor (T/Q) 2 . 

The situation is different in the case of = (r Poles of the Green's functions under integrals appear to 

be on the opposite sides of the real axis. As a result at uj — the integral diverges at finite x as widths 7~b and 
constants c tend to zero. Then one has to consider entire series to find T" 1 . After analysis of diagrams for and 
evaluation of their contribution to %p(w) we find that the most important diagrams within each order of g 2 are taken 
into account in the following equation: 



J-2 poo 

r i m J x + "> x ) = P mn + E ^^V^ / dyy\y\N(y)A(y)G l (x + y + u>)G* q (x + y) 

l,q J -°° 
! roo 

wS 7 mm S z nn / dyyN(y)A(y)r+- n (x + y + u,x + y)G m (x + y + u)G* n (x + y). (22) 



l,q 

h 2 T 2 

ttO 2 



—00 



Notice that this equation differs from that obtained in I within SWA. It takes into account only ladder diagrams 
shown in Fig. [5] whereas within SWA one has to consider also diagrams with crossing of two neighboring rungs. 
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Eq. 1221) can be easily solved if one notes that the last term in this equation should be taken into account only when 
W + c m ~ Cn\ <C T. Outside this interval the last term is much smaller than the second one. At such u the area of 
integration near the poles of the Green's functions is essential in the last term in Eq. (122[l and it gives the resonant 
contribution: 



Lq 



dyy\y\N{y)K{y)Gi{x + y + uj)G*Jx + y) 



2ih 2 S mm S* nn ( Z) (x + c n )N(-x - c n )A(x + c n ) - T /™ { " ■ % . (23) 



0/ " w + Cm - c n + ir 0m + ir an 

Eq. (|23[l has the same structure as the equation derived in I. The difference is that Po and the third term are of the 
order of g 2 rather than g 4 . Solving Eq. (|23|) one obtains 

T p mn ( x + LJ ^ x ) = p mn + —^^2S ml Pi q S% n d utl / dyy\y\N(y)A(y)Gi(x + y + u})G* q (x + y) 

l,q J -°° 

+ **'*^eA.e (i) 3 , +c ,^::u,„ ' p*) 

@ /T\ 3 

T m „ = ^(m-n) 2 - - . (25) 

We use in Eq. I|24|l that the following x will be important for calculation contributions to xp(uj) from T" 1 : |x+c n | <C T. 
Note, the third term in Eq. (|24l) is much greater than the second one when \u> + c m — c n \ <C io// 2 . It should be 
pointed out that F mn — for P — S z . Evidently, temperature dependences of r mn for P — S x ,S y and Jo given 
by Eq. (|15|l are the same. Within SWA the constant r mn has the same T-dependence but it is of the order of 

rkr ] oc / 4 T 3 /e 2 . 

IV. IMPURITY SUSCEPTIBILITY 

We can derive now the impurity susceptibility using the general expression J7J, Eqs. ©, I|1U|) and for the 
Green's function, Eqs. I|21|l and (|24|l for the branches of the vertex and Eq. for the number of pseudofcrmions 
N. It is convenient to discuss separately the cases of S = 1/2 and S > 1/2. 

A. Spin-1/2 impurity 

As a result of tedious calculations some details of which are presented in Appendix we have for the dynamical 
susceptibility of the impurity up to terms of the order of g 2 



1 2ir f 2 f°° , sgn(a;)A(x) 



PS^,P]d^ 

Rx 25T1 ' (2?) 



where we introduce the notation Y = Tr(Y), T = h 2 T 3 /(49 2 ), r = for P = S z and P = h 2 T 3 /{2& 2 ) for P = S' a: ^. 
Expression (|26|1 differs from that obtained in I only by the form of Po and P and by the thermal 1/s-correction to 
R x . The first term in Eq. I|26() is the Lorenz peak with the width P. The last term is the non-resonant part of the 
susceptibility. Its imaginary part at |cj| S> Po is proportional to sgn(w) and the real one contains the logarithmic 
singularity of the form ln(cij 2 + Pq). At T = and ui ^ the nonresonant contribution survives only and the 
susceptibility has the logarithmic singularity. It is clear that longitudinal host spins fluctuations become important 
if P > r^ swa) and P > These conditions are fulfilled when \g\/J <§C 5, i.e., for weakly coupled defects 

being considered. The contrast between results which do and do not take into account the longitudinal host spins 
fluctuations is illustrated by Fig.|SJ where we plot Imxj_(^)/w using Eq. 

The non-resonant term gives the main contribution to the susceptibility (|26|l when 

9 /P N 2 



lu ^> lu = — i — ) . (28) 
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Notice that lu oc f 2 T 2 within SWA. 

Using results of calculations presented in Appendix IE1 we find that static susceptibility Xp(0) does not depend on 
r and Jq: 



w = ^(l-4/ 2 UB x )+4R x; ^lii(|!) !«;, 



7T 3 /2 7 x |a;| tj-5/2 y (J2 _ j2)3/2 

It is taken into account in Eq. H3()[l that one cannot use Eq. (JHJ) for ImAi(s) because large x are important in the 
integral^ Eq. I|A2() is used to find U as it was done above for calculation of constants c. It is seen from Eq. I|29|) that 
the static susceptibility contains the free-spin- like term (4T) _1 which amplitude is slightly reduced by the interaction 
and the logarithmic correction proportional to f 2 . Static susceptibility has the same form as that obtained in I 
within SWA with the only difference in the form of the constant R x = R x swa) + [v{T)/(2s)]PS z [S z , P}. Notice that 
R x = R ( x swa) = 1/4 for P = S z and X||(0) = X^i )- In contrast R x = R ( x wa) + v{T)/(8s) = 1/8 + v(T)/(8s) for 
P = S x,y and the transverse static susceptibility has the form 



4T V 2 J J v^sB \sr)J 2tt9 \T 

Thus the longitudinal host spins fluctuations lead to another logarithmic term in x_l(0) 



B. Impurities with S > 1/2 

1. Dynamical susceptibility 

The calculations are slightly more tedious for impurity with S > 1/2. As a result of simple evaluations presented 
in Appendix iBl we lead to the following quite a cumbersome expression for the impurity susceptibility: 



^ ~t~ c m c n 4~ c 2i%xyi 



'I 2 (4n+4'i) ( £ 

2Re(P nm S n \ q P q iSj 1 n d nu ) f 2 
uj + c m - c n + ir 0n + ir 0m tt@ 



\p„^ i iK-:. - mr ) i^ + e^) r dx sgn(x)A( ^ — — 

/ 27t6 J_ oq x + uj + c m - c n + ir 0n + ir 0n 



dx\x\xN(x)A(x) 



e (c m +ir om )/T 



(x-u + ci- Cm ~ iFoi - ir 0m )(x + c q - c m + ir 0q - iPom) 
e (c n -tr on )/T 

(x + w + c q - Cn + ir Qq + ir 0n )(x + c; - c„ + iP 0n - ir i) 



(32) 



The first term here is zero for longitudinal susceptibility (P = S z and thus m = n). In the case of transverse 
susceptibility it gives 25 Lorenz peaks corresponding to transitions between impurity levels. Using Eq. (|17f) we obtain 
for the resonant frequencies 

4") = |cn - c n _i| = (2n - l)C(T)g 2 /J, (33) 

where n = S 1 , iS — 1, . . . ,n m ; n and n m i n = 1/2 and 1 for half-integer and integer spins, respectively. Then in contrast 
to half-integer spins there is no resonance peak at zero frequency in the case of integer spins. Notice that the distance 
between frequencies does not depend on S. The second and the third terms in Eq. (|32(l give contributions to the non- 
resonant term. The third term contains also a resonant contribution. As is explained in Appendix^] this contribution 
should be discarded because the resonant terms of the susceptibility are calculated within the zeroth order of g 2 . 

One makes sure after simple transformations that Eq. I|32|l coincides with Eq. (1201) if S = 1/2. It is more convenient 
to discuss expression (|32|) further in some limiting cases in which x( w ) nas simpler forms. Let us consider the 
susceptibility when temperature is greater than the impurity splitting, T 3> C(T)g 2 S 2 / J, and when T is much smaller 
than the energy between nearest impurity levels, T <C C(T)g 2 / J. 
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T C(T)g 2 S 2 1 J. Eq. 1M2I) can be simplified greatly in this case if one expands exponents e c "/ T up to the second 
term and notes that the nonresonant terms should be taken into account at |cj| 3> C{T)g 2 S/ J only. As a result one 
obtains 

/ \ 1 Vip 12 c m -c n +2ir nm f 2 f 00 sgn(x)A(x) 

T{2b + 1) ^— ' w + c m - c„ + 2ir„ m n®J-oo x + uj + id 

where i? x is given by Eq. I|27|). Notice that Eq. (|34[1 coincides with Eq. H26[) if 5 = 1/2. Imaginary part of the 
nonresonant term is greater than that of resonant terms when \u>\ 3> u)q, where ujq is given by Eq. I|28|l . The real part 
of the nonresonant term dominates in Eq. (|34[1 at |a>| ^> /O-i/G/T. To illustrate our results we plot in Fig. [7| the 
even-w function Imx±(w)/w using Eq. 134|) for S = 3/2 and S — 2. 

T <C C(T)g 2 /J. At such small temperatures only low- lying impurity levels contribute to susceptibility. The 
transverse susceptibility for integer impurity spin has the form 

Xl» = ~ W * 7 2|P01 f+ rr +4{\P i\ 2 (R^+R^)s g n(^c 1 -c ) 

uj + Co — c\ + ziIqi uj — (co — c\) + 2il oi B L V / 

- ""^ 

(1 + sgn(w + c g - c )) - ->■ -u;} > , (35) 



^ RefPomS^P^Sfjii^) (a; + c g — c ) 2 
^ uj + c m - c + ir 0m w + c 9 - ci + ir Qq + ir i 



where {uj — > — denotes terms inside the curly brackets with — u put instead uj. The first two terms in Eq. H35|) 
describe transitions between three low-lying impurity levels: {—1 <-* 0} and {0 <-> 1}. The last term in Eq. (|35() is the 
non-resonant part of the susceptibility that is zero between the resonance peaks uj — ±(co — ci). Beyond this interval 
it is negligibly small near the resonance peaks, i.e. at \oj + c — c\ | <C T 3 / 2 /\/0 and \uj — (c — ci)| -C T 3 / 2 /y/Q, but 
it gives the main contribution to imaginary part of the susceptibility outside these two intervals being of the order of 
f 2 . The real part of the non-resonant term is small compared to that of the resonant terms at \u>\ < 6. 
In the case of half-integer impurity spins the transverse susceptibility has the form 

, . . 1 1 2 2ir_ii \Pi 3 \ 2 \Pi 3 \ 2 

half—int/ \ 1 2 2 2 2 i 1 H 1 22 1 

X± ( LJ ) — 



T uj + 2ir_ii uj + ci — c 3 + 2iP\ s uj — (ci — C3 ) + 2iP\ . 

22 2 2 22 v 2 2 ' 2: 

•/ 2 ' 



n=±l 



X! I P »™! 2 + 4™) sgn(^ + c m - Cn) 



2Re{P nm S mq P ql Sld^) (uj + c q -c n ) 2 
- }^ — ; ^ ; V~7^ — r - ?^ 1 + s § n ( w + c <? ~ c ™)) ~ i w ~* _W A' ! Q) 

^ CJ + Cm ~ Cn + ll 0m + W On UJ + C q - Ci + tl o q + ll i 

m,q 7 L 

where the first three terms describe transitions between four low-lying impurity levels: {—3/2 <-> —1/2}, {—1/2 <-> 
1/2}, and {1/2 «-> 3/2}. The non-resonant part of the susceptibility is negligibly small in the vicinity of the resonance 
peaks, i.e. at \uj\ < T 2 /9, \uj + c 1/2 - c 3/2 | < T 3 / 2 /\/6 and |w - (c 1/2 - c 3/2 )| < T 3 / 2 /V&. Outside these three 
intervals the nonresonant term is of the order of f 2 and it gives the main contribution to imaginary part of the 
susceptibility. The real part of the nonresonant term is small compared to that of the resonant terms at \uj\ < O. 



2. Static susceptibility 

T ^S> C(T)g 2 S 2 1 J. Using results of Appendix 151 one finds for static susceptibility in this regime 



where R x and U are given by Eqs. (|27l) and (|30|l . respectively. Eq. (|37() coincides with Eq. I29|) at 5 = 1/2. It is seen 
that static susceptibility of large-5 impurities has the same structure as that of spin-^ defect. 

T <C C(T)g 2 1 J . At small temperature integer and half-integer impurity spins behave differently. Thus, transverse 
static susceptibility of an integer spin is a constant: 

Xl nt (0) - (38) 

Co - C\ 
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whereas longitudinal one is exponentially small: xj"*(0) « 0. 

Transverse static susceptibility of a half-integer spin has the form 



hal f—int (r.\ 

X± (0) 



1 

T 



P-Ll 

2 2 



2 \ f 2 /B 

- fUQ + 4Q^ In 



T 



Pi 3 
2 2 



Q — — I J]] | Pnm | 2 -Rj>j J]] PnmSmqPqlS^d^v J 
\n,m n,m,q,l I 



(39) 
(40) 



where primes over sums signify that summation is taken over projections 1/2 and —1/2 only. It is seen that to the 
first two terms in Eq. i|39|) contribute only transitions between states with projections ±1/2. The last term represents 
the contribution to the susceptibility from transitions between states with projections ±1/2 and ±3/2. It is of the 



order of l// 2 being much larger than the logarithmic term up to exponentially small temperatures T 
Longitudinal static susceptibility of a half- integer spin behaves like that of spin-i defect: 



Q e -i/r 



hal f - 




(41) 



OfJ. p Of J 
i J mq r ql J ln Ll . 



(42) 



where Q is given by Eq. Q4U|). In contrast to X± m *(0) there is no term proportional to l// 2 in x 
screening the logarithmic correction. 

As it is explained in I, there is a restriction on the range of validity of the resultant expressions for %p(o;). It is 
the consequence of the fact that the function ImA(w) has the form © if \u\ ^> r\ only. It is easy to see that in all 
calculations performed above one can use the function of the form © if the following condition on lu and widths Fq 
holds: max{min{F 0n }, \oj\} » r\. 



half—int 

•II 



(Q) 



V. INFLUENCE OF DEFECTS ON THE HOST SYSTEM 



We discuss in this section the influence of finite concentration of the defects n on the spin-wave spectrum of AF. 
It is demonstrated in I that in the vicinity of the points k = and k = ko, where ko is the antiferromagnetic vector, 
denominator of host spins Green's functions determined the spin-wave spectrum has the form 



X>(w,k) = u) 



2 



2 



^2 



(43) 



, \ 1 (kRioj 

«W = 2+ fe^' ( 44 ) 

where R12 is the vector connected host spins coupled to impurity (it is assumed for beginning that this vector is the 
same for all defects) and it is used that the unperturbed spectrum is linear at k <C ko and k ~ fco : £k = ck — \/8sJk. 
It is seen from Eqs. (|43|l and (|44|) that the spectrum depends on the direction of the momentum k as a result of 
interaction of magnons with the defects. This circumstance is a consequence of our assumption that the vector R12 
is the same for all impurities. In fact, it can have four directions and the value (Ri2k) 2 /fc 2 can have two different 
values: cos 2 </>k and sin 2 </>k, where </>k is the azimuthal angle of k. It easy to realize that it(k) = 1 if all four ways of 
coupling of the impurity with AF are equally possible. 

Let us turn to spin-i impurity. It is convenient to consider separately regimes of \ui\ ^> uio and \u\ <§; cjq, where u>o 
is given by Eq. I|28l) . In these cases the non-resonant and the resonant parts, respectively, dominate in the impurity 
susceptibility l|2f>(l . The results are summarized in Table [I] They are very similar to those derived in I within SWA. 
The only difference is in the form of wo, F and R x . It is seen that due to interaction with defects the spin- wave 
velocity acquires negative corrections and strong spin- wave damping arises proportional to nf 4 uj at u> ^> ujq. 

As it is explained in I the results obtained within our method have a limited range of validity which is also indicated 
in Table [I] This is because the interaction with defects modifies host spins Green's functions and the bare spectral 
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function ImA^°J (lv) given by Eq. renormalizes as follows: 



= ImA<°)( W )- S n/ 2 6Imxi( W )V, (45) 
B= i^/dkfi±^H)V«0.8. (46) 



7T 5 / 2 7 V J + </k 

The last term in Eq. (|45|) becomes larger than ImAjJ! (u;) at small enough energies indicated in Table [I] and the 
problem should be solved self-consistently. Corresponding consideration is out of the scope of the present paper. 

This large renormalization of the spectral function at small energies can be interpreted as a resonance scattering of 
spin waves on impurities. Similar situation exists in the case of asymmetrically coupled defects. Let us discuss, for 
example, AF with an impurity spin weakly coupled to one host spin. According to results obtained within T-matrix 
approach the defect leads to a resonant level inside the spin- wave bandii£*2ii2i The energy of this level is of the order of 
the energy cost for creation of spin excitation on the impurity. First perturbation corrections to observables increase 
rapidly near the resonant level due to the resonance scattering In the case of frustrated spin-i defect the energy of 
creation of the spin excitation on the impurity is zero and the resonance scattering takes place at zero energy. 

In the vicinity of the resonance peak, where our theory is not applicable, one might expect strong deviation from 
linearity of the real part of the spin-wave spectrum and an enhancement of the spin-wave damping as in the case of 
asymmetrically coupled defects^. 

Particular expressions for spin- wave velocity and damping in the case of S > 1/2 can be obtained straightforwardly 
using Eqs. QD, lj36*|) and |@3J. Because of their cumbersomeness we do not present here the results. We would like 
to point out only that there is the abnormal spin- wave damping proportional to nf (stemming from nonresonant 

part of XJ-M) when i) \uj\ > T 2 if T > C{T)g 2 S 2 /J, ii) \u>\ - w&l > T 3 / 2 for n ^ 1/2 and 



I, ,1 , 

CJ — U) res 



> T 2 if 



T <C C(T)g 2 / J, where are given by Eq. The range of validity of the results obtained within our approach 

can be found comparing the last term in Eq. (|45|l with the first one as it was done for S = 1/2. The results are 
not valid in the vicinity of the resonance peaks due to the resonance scattering of spin waves on impurities. At such 
energies one might expect strong deviation from linearity of the real part of the spin-wave spectrum and increase of 
the spin- wave damping. 

Thermodynamic quantities of AF with defects can be found, in principle, using impurity susceptibility. But due 
to the limiting range of validity of our results and due to the fact that thermodynamic quantities are expressed via 
integrals containing susceptibility one cannot calculate them because it is impossible to perform the integration over 
the essential energy area (in particular, over small energies). We explain this situation in detail in I by the example 
of the specific heat. We do not obtain any noticeable corrections to the specific heat after integration over energies 
at which our results are valid. 



VI. CONCLUSION 



In the present paper we discuss dynamical properties of an impurity spin coupled symmetrically to sublattices of 
2D Heisenberg quantum antiferromagnets (see Fig. ^) at T > (existence of a small interaction stabilizing the long 
range order at T ^ is implied). We continue our study on this subject started in I where spin-i defect was discussed 
and the host spins fluctuations are considered within the spin- wave approximation (SWA). In the present paper we 
i) go beyond SWA and ii) study impurities with spins S > 1/2. We show that the effect of the host system on the 
defect is completely described by the spectral function as it was within SWA. It is demonstrated that longitudinal 
host spins fluctuations lead to two important contributions to the spectral function that is proportional to oj 2 within 
SWA. They have the form r(T)T 2 ui / s and cv 2 v(T)/s, where s is the value of the host spin, r(T) and v(T) are series 
in (T/s) hi(T / [srj\) . First terms in these series are given by Eq. J5|. 

It is demonstrated that transverse spin-i impurity susceptibility X± i u ) nas the same structure as that obtained 
within SWA: the Lorenz peak and the non-resonant term. The difference is that the width of the peak r becomes 
larger being proportional to f 2 J(T/ J) 3 rather than f i J(T/ J) 3 . Besides some constants in the expression for x^^ ( w ) 
acquire thermal corrections. The width of the Lorenz peak in longitudinal impurity susceptibility is zero within our 
precision. The resultant expression for x^^H^) is given by Eq. (j2EJ). We observe that longitudinal host spins fluctua- 
tions lead to a new logarithmic correction to transverse static susceptibility l|31|) proportional to — {f 2 /[Js}) \n(T/[si]]). 
It is demonstrated that corrections to r from longitudinal host spins fluctuations dominate if \g\/J <C 5. 

We derive also expressions for x( w ) f° r symmetrically coupled impurities with S > 1/2. It is well-known that 
host spins fluctuations lead to an effective one-ion anisotropy on the impurity site in 3D AF: —C(T)S 2 g 2 / J, where 
C(T) > Oi 18 i 19 Then defects with S > 1/2 appear to be split and magnetization of sublattices is a hard axis for them. 



12 



We obtain the same one-ion anisotropy for defects in 2D AF. In particular we observe large reduction of the value 
of C(T) obtained within SWA by 1/s-corrections stemming from longitudinal host spins fluctuations. We show that 
thermal corrections reduce C(T) further and can change the sign of C(T) at T < TV. 

Then we find that if T is greater than the value of the defect splitting, transverse dynamical susceptibility contains 
2S Lorenz peaks corresponding to transitions between impurity levels and a nonresonant contribution. Widths of 
the peaks are proportional to / 2 J(T/J) 3 as in the case of S = 1/2. When T is lower than the distance between 
nearest impurity levels (~ C(T)g / J) there are only peaks corresponding to transitions between low- lying levels. 
Longitudinal susceptibility has only the nonresonant term. As it was in the case of spin-^ impurity, the imaginary 
part of the nonresonant term is a constant that leads to abnormal spin- wave damping proportional to nf 4 \ui\. At 
T > C(T)g 2 S 2 / J static susceptibility has the same structure as that for S = 1/2: term S(S + 1)/(3T) and the 
logarithmic correction. If T <C C(T)g 2 /J, x(0) behaves differently for integer (Eq. (j2HJ) an d half-integer (Eqs. 139|) 
and H41fl ) spins. 

Corrections to the spin- wave spectrum are obtained in the case of finite concentration of spin-i defects. They 
are summarized in Table [i] They have the same form as those obtained within SWA. The difference is that some 
constants acquire thermal corrections, T has another form and the quantity ujq separating two regimes appears to be 
proportional to T 2 /(sJ) rather than f 2 T 2 / J. In particular, we find strong spin-wave damping proportional to n/ 4 |a;| 
for all S originating from the non-resonant terms of the susceptibilities. 
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APPENDIX A: CALCULATION OF ImA M „(w) 

In this appendix properties of the imaginary part of the function A M „(a;) are discussed general expression for which 
is given by Eq. (J2J). We show that ImA^^) has the form JjJJ and find constants A and B, the characteristic energy 
O and tensor d M „. 

One has from Eq. (J2J 



A ""M = n £ [1 + cos ( kR i2)]<sVk}c (Al) 



k 

where (. . . ) u denote retarded Green's function, N is the number of spins in the lattice and R 12 is the vector connecting 
two host spins coupled to the defect. Thus we have to calculate the spin Green's functions (s^ k s%.) w . 

It is demonstrated in I that only diagonal components of (s_i c s k)w are nonzero. Only xx- and yy- components are 
nonzero within the spin- wave approximation (SWA). At T — we have for them in the isotropic 2D AF— 

ImA^) = -dua,^ J dk{[l + cos(kRi 2 )][J - Jk] + [l + cos([k + k ]Ri 2 )][Jo + Jk]} 

x-[(S( W - £k )-% + £k)], (A2) 
£k 

where d^ v — 5^(1 — 5p Z ), ko is the antiferromagnetic vector, the lattice constant is taken to be equal to unity, 
Jk = 2J(cosk x + cosfcj,), ek = s yj Jq — is the spin- wave energy and the integral is over the magnetic Brillouin 
zone. If <C sJ we have Jk ~ Jo — Jk 2 , cos(kRi2) » 1 — (kRi 2 ) 2 /2 and ek = ck = y/8sJk. Notice that 
(koR.12) = Ti" mod 27T if the impurity is coupled to spins from different sublattices and both terms in the curly 
brackets in Eq. i|A2|) are proportional to k 2 . Then integration in Eq. (|A2|) can be easily carried out if one takes 



advantage of the approximation for magnons similar to Debye one for phonons: the spectrum is assumed to be linear, 
£k = ck, up to cut-off momentum fee defined from the equation 4irN — V \ Q e dkk, where V is the area of the lattice. 
As a result we lead to terms proportional to d^ u {l — <5 MZ ) in expression J3J for ImA M „(cL>) in which 

6 = ck e = sJSVtt, (A3) 
A = % (A4) 
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The constant A should be multiplied by 2 if the defect is coupled to four host spins (two by two from each sublattice). 

It is easy to conclude that if a small interplane interaction of the value of r\ <C J is taken into account the above 
result for ImA^^u) is valid when \u>\ ^ 77. At the same time ImA^fw) has another w-dependence if \u>\ < r\. It is 
well established2^Si2L22, that spin waves are well defined in paramagnetic phase of 2D AF if their wavelength is much 
smaller than the correlation length £ oc exp(const/T). Thus, the above result for ImA A1 y(cj) is valid when |u;| Ja/£. 

Let us go beyond SWA. It is well known that the spin-wave interaction gives negligibly small corrections to the 
spin-wave spectrum and other physical quantities in isotropic Heisenberg J^Ysm^^^l^^ For instance, simple 
calculations give that transverse components of ImA M „(u;) given by Eq. (|A2|) are multiplied by 




1 \ -> £k 
1 > — 

k 



0.08 /A , 

1 - — (A5) 



after taking into account first 1/s-corrections. Meanwhile small interactions of the value of r\ possibly could lead to 
a great renormalization of physical observable quantities at <C rj^^ At the same time we derive ImA^^fw) at 
3> r/. Thus we can use the expressions for xx- and yy- components obtained above within SWA. 
We have to take into account also the longitudinal spin susceptibility (s k sl k )w ■ The first correction to it is of the first 
order of 1/s. This quantity has been examined in Ref^. Nevertheless we derive now the corresponding expressions 
in a more convenient for us form. Using retarded Greens functions <j(w,k) = (akjS^, /(w,k) = (a k ,a_k) w , 

g(cu, k) = (a^ k ,a_k)w = <?*(— — k) and /^(w,k) = (a^L k ,a k ) w = /*(— k>, — k), where ak and a k are operators of 
magnons creation and annihilation, we have 

(s k s -kL = -— ^2 X] [/(«wi,ki)/ t (iw-iwi,k2)+g(iwi,ki)5(iw-iwi,k 2 )]. (A6) 

0J1 ki+k 2 =k-|-ko 

Functions / and g were calculated, e.g., in I. After simple evaluations one obtains 
7r <r-^ 1 — cos([ki + kajRia) 



2N ' 2 k k ekl£k2 

ki ,k2 



ImA zz (tij) 

ki,k2 

x |[1 + 2iV(e kl )] [(sJ ) 2 + s 2 J kl J k2 - e kl e k2 ] [5(u> + e kl + e k2 ) - S(uj - e kl - e k2 )] 
+ 2[AT(e k3 ) - AT(e kl )] [(s J f + s 2 J kl J k2 + e kl £k 2 ] S{uj + e kl -e k2 )|, (A7) 

where the summations are over chemical Brillouin zone. The first term in the curly brackets in Eq. i|A7(l corresponds 
to emission or absorption of two magnons whereas the last term describes scattering of one spin wave. Let us discuss 
firstly the case of T = 0. The Planck's functions in Eq. (|A7I) are zero and ImA zz (w) is proportional to cu 3 . Therefore, 
it is small in comparison with xx- and yy- components obtained above. 

At finite temperatures we have two important contributions to ImA zz (w). One of them arises at \u>\ T and 
originates from the second term in the curly brackets in Eq. (|A7|) . It can be brought to the form — ujT 2 B /(s<d 3 ), 
where 

„ 240F r , x 2 e x 8vr 5 / 2 

The second contribution originates from both terms in the curly brackets in Eq. (|A7I) and has the form 
— Aw 2 sgn(w) , u(T)/(s9 2 ), where v(T) — 8sJ/iV^ k A^(e k )/ek. The logarithmic infra-red singularity in this expres- 
sion is screened by the interaction stabilizing the long range order at T ^ 0: v(T) = (8/v / 7r)(T/0) ln(T / [sr/]) . There 
is no parameter screening this singularity in isotropic 2D AF. This singularity signifies that spin- wave approach does 
not suite for discussing the impurity dynamics in isotropic 2D AF. 

It should be stressed that higher order 1/s-corrections can give contributions proportional to products of ui 2 (uT 2 ) 
and powers of [T/ (sQ)] ln(T/[srj\). As a result the obtained above correction proportional to loT 2 should be multiplied 
by a function r(T). This function as well as v(T) appear to be series in [T/ (sQ)] ln(T/[sr)]). To find general expression 
for r(T) and v(T) is out of the scope of the present paper. Thus, we lead to zz-component of ImA^fw) in Eq. l|3]l. 

Notice that in the case of impurity coupled to one host spin, terms in Eq. (|A2(I containing cos(kRi2) should 
be discarded and the spectral function appears to be proportional to a constant within SWA. There is no large 
corrections to this constant from higher order 1/s-terms. Hence, one could not expect that longitudinal fluctuations 
play any remarkable role in this case. Then, we confirm that the so-called T-matrix approachii&ii is appropriate for 
asymmetrically coupled defects. Recall that it bases on Dyson-Maleev transformation and deals with only bilinear 
part of the Hamiltonian. 
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APPENDIX B: CALCULATION OF THE IMPURITY SUSCEPTIBILITY 



We present in this appendix some details of the impurity dynamical susceptibility calculation. We use for this 
general expression (JJJ and Eqs. (JSJ, (|21[1 and l|24|) for the Green's function and the branches of the vertex. The 
expression for Xp( w ) i s derived up to the order of g 2 . 

According to Eqs. J7J, (|21(l and (|24fl the dynamical susceptibility can be represented as a sum of three components. 
The first one, xi( w )j originates from Eq. J7|) as a result of replacement of the vertex by unity. The second, X2(k-0> 
appears from / 2 -terms in Eqs. (|21(l and (|24(l . The third, X3( w )> is a result of replacement of the vertex by the third 
term from Eq. H24|) . Calculations of these three quantities are similar to those performed in I. Then we discuss the 
results only. 

For xi( tL ') we obtain 



xiM 

e-Vr 
Af 
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En 
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(1) 



r- 
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(Bl) 



For X2(w) one has 



X2M 
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Eq. (B2j has a simple form at |w| > C{T)g 2 S/J: 



01 > 



X2H = - 

F° r X3( w ) one obtains 
e -A/T 



-A/T 
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X3M = 
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(B2) 



(B3) 



(B4) 



It should be stressed that we use while calculating Eq. (|B4|> that lo is close to one of the resonance frequencies lo. 
given by Eq. 



LO 



.,(") 
LUres 



<C T. Meanwhile we can use expression ()B4|) at all lo as it is much smaller than x\{oo) 

far from the resonances. 

It should be noted that the resonant terms in xi, 2(^-0 and X3( w ) are calculated in the order of g 2 and g°, respectively. 
We would like to stress that calculation of X3( w ) m higher orders demands taking into account in equation l|22|) for 
^~Pmn( x + cc) not only the most singular diagrams in each order of g 2 . Their analysis is a cumbersome task that is 
out of the scope of the present paper. Thus we restrict ourself in this paper by calculation of resonant terms in the 
dynamical susceptibility in the order of g° and discard the resonant terms in Xi( w ) an d X2O- 1 -') that are of the order 
of g 2 . Notice that one can keep these terms while calculating static susceptibility because Xa(0) = an d the problem 
of the cumbersome analysis of the higher order diagrams does not arise^ 
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TABLE I: Renormalization of spin-wave damping and velocity by interaction with spin-j impurities in two regimes: 2> ljo 
and \u\ <C u>o, where luq is given by Eq. 1281 . Ranges of validity of the theory in each regime are also presented (see the text). 





M > 


\ui\ -C <^o 


spin-wave damping 
spin-wave velocity 
range of validity of the theory 




e^!r^n£l u(k) 
T uj + 4r 87r 


c^l-ln 1© ^^(k) 

y » o.iv 'R^V^f 


l^^ 2 + 4r2 )»o.oom/^ 




FIG. 1: Unit cell of 2D AF with impurity spin coupled (a) symmetrically and (b) asymmetrically to AF sublattices. Strengths 
of coupling with corresponding host spins g and <?i 7^ <?2 are depicted. The local Neel order is also shown. Only symmetrically 
coupled impurities are discussed in the present paper. 
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FIG. 3: Lower-order diagrams for pseudofermion self-energy part. 



a) 



b) 




FIG. 4: Examples of diagrams for pseudofermion self-energy part (a) and vertex (b) containing vertexes with two and three 
magnons lines. Such diagrams are negligibly small compared to those presented in Figs. |3]and 



Yp(x + LJj x) = 



> + > + > + 3^ 
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FIG. 5: Lower-order diagrams for pseudofermion vertex Tp(x + u, x) 



+ . 



FIG. 2: Lower-order diagrams for the impurity dynamical susceptibility \p(uj) and a graphical representation of the result 
of the overall series summation. Lines with arrows and wavy lines represent pseudofermion and magnons Green's functions, 
respectively. 
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co/@ 

FIG. 6: The even-a; function lm(x±(w)/w)0 2 for spin-^ impurity at T/Q = 0.1 and / = 0.1 with longitudinal host spins 
fluctuations being (solid line) and not being (dotted line) taken into account. 
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FIG. 7: The even-w function lm(x±(w)/w)0 2 calculated using Eq. (1341 for impurity with S = 3/2 (solid line) and S = 2 
(dotted line) at T/Q = 0.05 and / = 0.1. 



